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An extensive computational effort has been performed in order to investigate the nature 
of unsteady flow in the fuel line supplying the three Space Shuttle Main Engines during 
flight Evidence of high cycle fatigue (HCF) in the flow liner one diameter upstream of the 
Low Pressure Fuel Pump inducer has been observed in several locations. The analysis 
presented in this report has the objective of determining the driving mechanisms inducing 
HCF and the associated fluid flow phenomena. The simulations have been performed using 
two different computational codes, the NASA MSFC PHANTOM code and the Pratt and 
Whitney Rocketdyne ENIGMA code. The fuel flow through the flow liner and the pump 
inducer have been modeled in full three-dimensional geometry, and the results of the 
computations compared with test data taken during hot fire tests at NASA Stennis Space 
Center, and cold-flow water flow test data obtained at NASA MSFC. The numerical results 
indicate that unsteady pressure fluctuations at specific frequencies develop in the duct at the 
flow-liner location. Detailed frequency analysis of the flow disturbances is presented. The 
unsteadiness is believed to be an important source for fluctuating pressures generating high 
cycle fatigue. 


Nomenclature 

N = rotational speed 

NSS - suction specific speed 

Pi - inlet static pressure 

Q ” volumetric flow rate 

RPL = rated power level 

RPM = revolutions per minute 

L Introduction 

D uring a routine post-flight inspection of the liquid hydrogen (LH2) feed lines in the Shuttle Main Propulsion 
System (SMPS) leading to the Space Shuttle Main Engines (SSME), cracks were discovered between adjacent 
slots in a flow liner just upstream of the low pressure fuel pump (LPFP) inducer (see Fig. 1). The flow liner is 
located at a gimbal joint less than one duct diameter upstream of the LPFP inducer and has two rows of 38 slots each 
(see Fig. 2). The slots are stamped into the liner to allow for cleaning of the bellows located in the cavity behind the 
liner. The majority of the cracks were found in the downstream slots (closest to the inducer). Upon further 
examination it has been found that the cracks propagated in both the flow and slot-to-slot directions. 

An extensive series of hot-fire tests, water-flow tests, air-flow tests and numerical simulations were performed to 
determine the cause of the cracks, as well as to determine the risks they posed. The purpose of this paper is to 
present results from some of the numerical simulations performed for the feed line and LPFP inducer. The numerical 
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simulations were conducted using the two different analyses, the PHANTOM code of NASA Marshall Space Flight 
Center [1] and the ENIGMA code of Pratt & Whitney Rocketdyne [2]. The numerical simulations were performed 
for LH2, water (H20) and air. Both engine (hot fire) and scaled (water, air flow) operating conditions were 
considered. The effects of eccentric tip clearance and duct angulation were also studied. The general goal of the 
simulations was to characterize the unsteadiness upstream of the SSME LPFP inducer, while the specific goals were 
to isolate the sources of unsteadiness reaching the flow liner and provide guidance for the tests. The predicted results 
were compared with analytical solutions, hot-fire data and water flow data. 



Figure 1. Slot-to-slot and flow stream direction cracks in the flow liner. 



Figure 2. Flow liner at gimbal joint, just upstream of the LPFP inducer. 
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II. Numerical Procedures 

The governing equations in the PHANTOM code are the three-dimensional, unsteady, Navier-Stokes equations. 
The equations have been written in the Generalized Equation Set (GES) format, enabling it to be used for both 
liquids and gases at operating conditions ranging from incompressible to supersonic flow. The code employs a 
system of overset O-grids and H-grids to discretize the flow field. The grids move to simulate blade motion. A 
modified Baldwin-Lomax turbulence model is used for turbulence closure [3]. The code contains two options for the 
fluid properties. The first option is based on the equations of state, thermodynamic departure functions and 
corresponding state principles constructed by Oefelein [4]. The second option is based on splines generated from the 
NIST Tables [5]. Message Passing Interface (MPI) protocols are used for parallel simulations. The code can model 
two-phase flows (e.g., cavitation) using the GES extensions of Merkle et. al. [6] and Ahuja et al. [7]. A detailed 
description of the code/algorithm development, as well as its application to several turbine and pump test cases, is 
presented in Refs. 1 and 8. 

The Enigma CFD code is a structured-grid finite-difference code capable of solving incompressible and subsonic 
compressible, steady (e.g., frozen-rotor) and unsteady (e.g., rotor-stator) flow problems and it features a parallel 
solution algorithm implemented in terms of the MPI parallel programming library. The numerical method uses the 
standard k-s turbulence model with wall functions and also features a length-scale filtered k-s. turbulence model 
suitable for capturing resident flow unsteadiness. More details of the code and its development can be found in Refs. 
2, 9 and 10. 


III. Numerical Simulations 

The SSME LPFP inducer operates in LH2 and has four main blades and four splitter blades. At the nominal 
(104.5% RPL) operating conditions the inducer rotates at 15,761 RPM and has a mass flow of approximately 155 
lbm/sec. Simulations were performed for operating conditions from 65% RPL to 104.5% RPL. 

A. Phantom Simulations 

Several sets of simulations were performed with the PHANTOM code. The simulations were performed at 
engine scale conditions for LH2, H20 and air over the range of flow conditions between 65% and 104.5% RPL. A 
series of CFD simulations for H20 were also performed at 51% scale, corresponding to the actual hardware used in 
the water flow tests [11]. The full-annulus computational grid for the inducer (including 22 inches of the feed line) 
contained approximately 2.8 million grid points for cases without cavitation modeled, and 4.6 million grid points for 
cases including cavitation. The grids for the non-cavitating simulations are shown in Fig. 3. It was decided that only 
the main blades would be used in the simulations because the splitter blades are located towards the downstream end 
of the inducer passage, and would have little effect on the unsteadiness at the flow liner. In addition, the inducer 
blades were allowed to scrape along the shroud (i.e., the tip clearance was set to 0) in the non-cavitating simulations. 
In the simulations with cavitation the nominal value for the tip clearance was used. Note that the assumptions made 
on the splitter blades and tip clearance (in the non-cavitating cases) will have an impact on the performance of the 
inducer. The simulations were performed on 12 processors of the SGI Altix clusters located at NASA Ames 
Research Center and NASA Marshall Space Flight Center. Approximately 6 revolutions of the inducer equate to one 
through-flow cycle (i.e., the time required by a fluid particle to move from the computational inlet to the 
computational exit), so all the non-cavitating simulations were run for at least 12 revolutions. The cavitating 
simulations were initiated from the non-cavitating solutions and run for an additional 5 revolutions. 
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Figure 3. Computational grid for the non-cavitating PHANTOM simulations. 


Two views of the time-averaged static pressure contours from the 104.5% RPL simulation in LH2 are shown in 
Figs. 4 and 5. The contours highlight the low pressure region near the tip leading edge where the inception of 
cavitation is normally observed. The predicted time-averaged static and total head rise coefficients are compared 
with results of half-scale water flow tests in Fig. 6 over a range of operating conditions [11]. In general, fair 
agreement is noted between the predicted and experimental performance curves. The differences are attributed to 
several factors, including: a) the baseline simulations were assumed to be non-cavitating, and b) the baseline 
simulations were performed without the inclusion of tip clearance 



Figure 4. Static pressure (psi) on the SSME LPFP inducer - LH2 - view 1. 
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Figure 5. Static pressure (psi) on the SSME LPFP inducer - LH2 - view 2. 

The large difference in static pressure between the suction and pressure surfaces (i.e., loading) at the inducer tip 
leading edge creates a movement of flow upstream from this region. This “backflow” from the inducer is composed 
of two distinct phenomena. The first is a convective component, which behaves much like the wake from an airfoil. 
It is characterized by swirling flow whose axial flow component moves upstream towards the inlet. The second 
component is a potential or pressure wave. The pressure wave component, in some cases, can persist much longer 
than the convective component. Figure 7 shows contours of the axial velocity (ft/sec), which highlights the 
backflow, for 104.5% RPL conditions. Figure 8 displays the trends in the backflow for LH2 and H20 over a range 
of operating conditions. The convective component of the backflow from inducer reaches the flow liner (7 inches 
upstream of the inducer) even at 104.5% RPL, and encompasses the flow liner near 80% RPL. 



Figure 6. Predicted and experimental head rise coefficients for the LPFP inducer. 
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Figure 7. Contours of axial velocity highlighting backflow from the LPFP inducer 

The predicted attenuation of the unsteady pressure waves upstream of the inducer at 65% and 104.5% RPL, 
along with the results of an analytical wave model [12] and hot-fire data [13] are shown in Figs. 9 and 10. Also 
shown in Figs. 9 and 10 is a decomposition of the unsteady pressure at various locations upstream of the inducer 
leading edge (X=0 in) for 104.5% RPL. The unsteady pressure amplitude at the inducer leading edge predicted in 
the CFD was used as input to the wave model. All four data sets show an exponential decrease in the pressure with 
distance upstream of the inducer leading edge. The rate of decay is somewhat greater in the CFD prediction than the 
data indicates. 

The amplitude of the unsteady pressure is low, but not insignificant, by the time it reaches the flow liner. The 
predominant unsteadiness occurs at the inducer blade passing frequency (approximately 1000 Hz) or its first 
harmonic. Many structural frequencies lie in the neighborhood of 1000 Hz and may contribute to the cracking 
process. The energy content predicted by the CFD simulations matches the energy content measured in the tests, 
however, the test data spectrum displays a significant amount of additional frequencies in terms of higher harmonics 
and noise. The additional frequencies identified in the experiments are due to the presence of phenomena such as 
cavitation, which were not modeled in the simulations shown in Figs. 9 and 10. 



Figure 8. Variation of backflow with power level. 
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Figure 9. Decay of the unsteady pressure upstream of inducer (X=0 in is the inducer leading edge) - LH2 


CFD simulations of water flow LPFP and duct operation have been made for full-scale hardware and half-scale 
hardware, the latter case corresponding to the actual water flow test hardware. In each case, a flow scaling has been 
made in order to obtain the flow regime equivalent to 104.5% and 65% RPL. In the case of full scale hardware, an 
arbitrary rotational speed was chosen, and the flow coefficient value matched. The predicted unsteady pressure data 
for the full scale water simulations at 104.5% RPL is compared to half-scale water rig data in Figs. 1 1 to 13 [13]. 
The CFD predicted decay rates show good agreement with the measurements. 

The CFD simulations performed for water at the exact scale of the hardware used in testing allows for a more 
consistent comparison between prediction and measurement. In the water flow experiments, the flow regime 
corresponding to the 104.5% RPL of the SSME has been considered, i.e., the 100% Q/N regime [1 1]. Consequently, 
the 65% RPL of the SSME corresponds to 50% Q/N regime in water flow testing. Figure 12 presents a direct 
comparison between CFD and test in terms of pressure fluctuations upstream of the inducer for 100% Q/N, while 
Fig. 13 shows the comparison for 50% Q/N. For 100% Q/N, the main frequencies identified in both the simulation 
and the experiment include: the blade passing (4N) and 2 upper harmonics (8N and 16N). The predicted amplitudes 
show extremely close agreement to the experiment. The pressure wave attenuation is achieved at half a diameter 
upstream, similar to the Fig. 11. For 50% Q/N the CFD prediction displays the 4N and 2 upper harmonics, while the 
measurement indicates a full train of upper harmonics frequencies. As mentioned above, cavitation processes are a 
possible cause for the richer spectrum, and stronger amplitudes in the experiment. The decay rates for the pressure 
fluctuation upstream is nevertheless well predicted: it is remarkable that in spite of more than double the amplitude 
at inducer leading edge, the pressure wave depth of penetration upstream is not larger for the 50% Q/N case. Both 
measurement and CFD prediction show this clearly. For both regimes, a low frequency corresponding to IN 
propagates upstream with less attenuation. Both the numerical and experimental data indicate about 0.5 psi 
amplitude for the IN wave. 
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Figure 10. Decomposition of the unsteady pressure upstream of the inducer - LH2. 



Figure 11. Comparison between unsteady pressure predicted for the full-scale rig with the experimental 

data from a half-scale rig - H20. 
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Figure 12. Comparisons between unsteady pressure predicted in half-scale H20 simulation and half-scale 
rig data. Test data FFT is marked in red at 4 locations - 100% Q/N. 



Frequency/N 


Upstream distance UD 


Figure 13. Comparisons between unsteady pressure predicted in half-scale H20 simulation and half-scale 
rig data. Test data FFT is marked in red at 4 locations - 50% Q/N. 

One of the features of the feed line is that it can angulate at the gimbal joint as the engine is rotated in flight. 
Simulations were performed with the inducer angulated 3 degrees with respect to the feed line. The value of 3 
degrees was based on estimates of the maximum deflection during flight. The objective of the simulations was to 
determine if the angulation generates a spatial 50 disturbance which could excite the flow liner structure. Figure 14 
and 1 5 contain contours of the time-averaged static pressure (psi) and velocity (ft/sec). The angulation in the duct is 
clearly visible in this figure. 

The angulation in the duct causes the resulting pressure and velocity fields to be asymmetric upstream of the 
inducer. The asymmetric backflow results in the flow liner experiencing different amounts of unsteadiness 
depending on the circumferential location. Figure 16 shows spatial decompositions of the pressure at the top and 
bottom sides of the duct near the angulation location. The resultant spatial frequency at 50 predicted in the CFD was 
small, and this fact was later verified in water flow tests [12]. Simulations were also performed at 65% and 
80%RPL, as well as in H20. All the simulations showed small levels of unsteadiness at 50. 
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Figure 14. Time-averaged static pressure (psi) for 3-degree angulation simulation - LH2 - 104.5% RPL. 



Figure 15. Time-averaged axial velocity (ft/sec) for 3-degree angulation simulation - LH2 - 104.5% RPL. 

Simulations were also performed to study the onset and growth of cavitating flow regions as the inlet pressure 
was decreased. The simulations were performed for H20 and were designed to model the experiments of Skelley at 
al. [11]. The computational grids were refined in the spanwise and streamwise directions for the cavitation 
simulations, for a total of 4.6 million grid points. In addition, the nominal tip clearance was added to the model. 
Figures 16 and 17 contain plots of the 0.20 vapor fraction iso-surface for 4 values of the suction specific speed. The 
onset of the cavitation process is characterized by small regions of cavitation being entrained in the backflow 
originating from the leading edge tip region of the inducer blades. Increasing the suction specific speed causes the 
cavitation region to extend further upstream of, and in the inter-blade region downstream of, the inducer leading 
edge. Further increasing the suction specific speed causes the cavitating flow from the individual blades to coalesce 
and form a cylinder of cavitation upstream of the leading edge. In addition, the cavitation penetrates further radially 
into the core flow. Figure 18 shows a comparison of the predicted and experimental values of the total head 
coefficient as a function of the suction specific speed. The bulk differences in the head coefficient values are caused 
by the omission of the splitter blades in the simulations. The predictions show the head coefficient fall-off at a lower 
suction specific speed than indicated by the experimental data. The more rapid fall-off is probably also caused by 
omitting the splitter blades in the simulations, as well as the need to tune some of the parameters in the cavitation 
model. Additional simulations are being performed to more accurately determine the breakdown point in the CFD. 


10 

American Institute of Aeronautics and Astronautics 



Figure 16. Vapor fraction of 0.20 iso-surface for the LPFP inducer: a) Pi=45 psi, NSS=5025, b) Pi=15.2 psi, 

NSS=1 1,500. 



Figure 17. Vapor fraction of 0.20 iso-surface for the LPFP inducer: a) Pi=8.7 psi, NSS=17,775, b) Pi=4.1 

psi, NSS=32,725. 


B. ENIGMA Simulations 

Unsteady flow simulations were performed to examine the effect of an elliptical housing and eccentric housing 
on the pressure perturbations propagated to the downstream flowliner slots. The rationale for looking at these 
configurations was that it was hypothesized that a non-uniform tip clearance would generate high-order spatial 
modes that could interact with the flowliner slots. The flow model includes a suction inlet duct and housing and the 
4-full blade/4-partial blade inducer with non-uniform tip clearance. Figure 19 shows a schematic of the flow 
domain. The flow grids were created using the APPT grid tool [15]. The flow model contains about 1 .2 million flow 
nodes and includes the duct and the full-360 inducer. 

The flowliner slots were not modeled but enough of the suction inlet duct was modeled to envelop the location of 
the slots. Figure 20 shows the overlapping grid system used and the location of key geometrical elements. 
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Figure 18. Predicted and experimental values of total head coefficient as a function of suction specific speed. 


The elliptical tip clearance model has a spatial period of 1 80 degrees with a min/max tip clearance distribution of 
0.004”/0.0r\ The elliptical tip clearance represents a centered rotor within a housing that is not circular. The 
eccentric tip clearance model has a spatial period of 360 degrees and a min/max tip clearance distribution of 
0.004”/0.0r\ The eccentric clearance model represents the situation when the rotor is not centered with the 
housing. 

Unsteady results were calculated for 104.5% RPL conditions with LH2. The calculations utilized a sliding mesh 
with the inducer blade hub and blade mesh rotating relative to the fixed duct and housing. The numerical time-step 
was selected so that the numerical Nyquist frequency was about 23641 Hz. Computed pressure time signals and 
fast-Fourier transforms (FFTs) are shown in Figs. 21 and 22. Figures 21 and 22 also show the generation of the 
blade passage frequency (N) harmonics, 4N, 8N, 12N etc. The character of the harmonics is similar for both the 
elliptical and eccentric tip clearance cases. 

The interaction of inducer blade perturbations with the non-uniform tip clearance is illustrated in Figs. 23 and 24. 
The spectrum shows the spatial (0) modal content, at an instant of time, of the circumferential pressure distribution 
near the downstream slots. The pressure distribution contains about 100 points. The character of the harmonics is 
different for both the elliptical and eccentric tip clearance cases. Because of the symmetry of the elliptical tip 
clearance distribution, a 20 mode is generated due to the tip clearance-blade interaction (Fig. 23). For the eccentric 
tip clearance simulation, a 10 mode and a very weak 50 are generated due to the tip clearance-blade interaction (Fig. 
24). 


IV. Conclusions 

A series of numerical simulations have been performed for the feed line and low pressure fuel pump inducer of 
the Space Shuttle Main Engine in an effort to determine possible fluid dynamic mechanisms generating cracks in a 
flow liner upstream of the inducer. The results of a series of the numerical simulations in liquid hydrogen and water 
showed close agreement with results of hot-fire and water flow tests. The simulations and experiments indicate that 
unsteadiness generated by the low-pressure fuel pump inducer probably plays a major role in the generation of the 
cracks, as the dominant frequency of the unsteadiness is in the vicinity of several structural frequencies. Angulation 
of the duct upstream of the inducer and asymmetric tip clearance do not appear to contribute directly to the crack 
generation process. 
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Figure 19. Perspective view of flow model for elliptical and eccentric tip clearance study. 



Figure 20. Duct and inducer meridional grid plane: d/s indicates the location of the downstream slots relative 

to the flow model. 
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Figure 21. Elliptical tip-clearance distribution: computed pressure time signals and spectrums near tip 



Figure 22. Eccentric tip-clearance distribution: computed pressure time signals and spectrums near tip 
leading edge and near flow liner downstream (D/S) slot location. 
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Figure 23. Generation of 20 with elliptical tip clearance model. 
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Figure 24. Generation of weak 50 with eccentric tip clearance model. 
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Numerical Methods 



• Two codes used to perform the simulations 

•Phantom (MSFC) 

•Enigma (PW - Rocketdyne) 

• Phantom 

• 3D unsteady finite-difference code 

• Generalized Equation Set 

• liquids, gases and two-phase flows 

• incompressible thru supersonic flow 

• Moving grids 

• Enigma 

• 3D unsteady finite-difference code 

• Incompressible thru compressible subsonic flow 

• Moving grids 

• Codes produce similar results on similar test cases 


• Geometry 

• SSME low-pressure fuel pump inducer 

• includes ducting upstream of inducer 

• 4 main blades and 4 splitter blades 

• 65% to 104.5% RPL 

• Configurations 

• Air rig (MSFC) - to investigate acoustics, 51% scale 

• Water rig (MSFC) - steady/unsteady pressure, 51% scale 

• Hot-fire test (SSC) - LH2, engine hardware 

• Effects considered 

• power level 

• inducer tip clearance eccentricity 

• angulation of duct 

• operating fluid (LH2, H20, air) and two-phase flow 
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Numerical Simulations 
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Phantom Simulations - 1 




• Full 360 degree simulation 

• splitter blades neglected 

• 22 inches of upstream duct modeled 

• full-scale (hot-fire) and 51% scale (H20, air) models 

• 2.8 million and 4.6 million grid points 
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Phantom Simulations - II 


Time-averaged static pressure (psi) for LH2 at 104.5 RPL 


p24G 000 
® 0 000 


• Large loading near the tip leading edge -> leads to backflow 

• Low pressures noted on suction side of the tip leading edge -> location where 
inception of cavitation usually occurs 
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Phantom Simulations - III 




Comparison of predicted head coefficients in air, H20 and LH2 
(all full-scale) with data obtained in the 51% scale water-rig tests 
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Phantom Simulations - IV 


Trends in backflow upstream of the inducer 




Axial velocity (ft/sec) for 104.5% RPL - LH2 



Backflow from inducer reaches the downstream slot for all flow conditions, 
and envelopes the entire flow liner between 65% and 80% RPL 
(note: more cracks were observed in the downstream slots) 
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Phantom Simulations - V 



Decay of the unsteady pressure waves upstream of the inducer - LH2 



• The predicted decay matches well with the results of an analytical wave model and hot-fire data 

• The unsteadiness is reduced, but not insignificant, when is reaches the flow liner 

• The predominant unsteadiness is at the inducer blade passing frequency (or a harmonic) 
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Phantom Simulations - VI 


Unsteady pressure decay predicted in full-scale H20 simulation 
compared with data from 51% scale water rig tests 




•Very good agreement between predicted results and test data 
• Predominant frequency at inducer blade passing frequency or harmonics 
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Phantom Simulations - VII 


Unsteady pressure decay predicted in 51% scale H20 simulation 
compared with data from 51% scale water rig tests 



Frequency^* Frequency/N 


100% Q/N 50% Q/N 

• Predicted frequencies and amplitude show very good agreement with water rig data 

• Unsteady pressure waves reach the flow liner 
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Phantom Simulations - VIII 


Effects of 3 deg duct angulation (engine movement) 



Time-averaged static pressure (psi) Time-averaged axial velocity (ft/sec) 




• 3 degree angulation corresponds to maximum value expected during flight 

• Duct angulates at approximately the midpoint of the flow liner 

• Concern is the generation of unsteadiness at 50 

• The backflow is asymmetric on the top and bottom sides of the duct 

• The upstream propagating pressure waves are asymmetric on the top 
and bottom sides of the duct 
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Phantom Simulations - IX 


Spatial unsteadiness at two instants in time 




• Angulation does not induce significant 56 unsteadiness (confirmed in water flow tests) 
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Phantom Simulations - X 














Vapor fraction increases as the inlet pressure decreases -> 


# 






Nss=5025 Nss=1 1,505 Nss=1 7,775 Nss=32,725 
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Phantom Simulations - XII 




Suction Specific Speed, NSS 
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Enigma Simulations - 1 



• Full 360 degree simulation 

• splitter blades included 

• ~25 inches of upstream duct modeled 

• full-scale model, eccentricity effects included 

• 1.2 million grid points 

• LH2 
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Enigma Simulations - II 



Elliptical tip-clearance distribution 
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Enigma Simulations - III 



Elliptical tip-clearance distribution 
Near D/S Slots 




Blade/clearance interaction generates 20 unsteadiness 
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Enigma Simulations - IV 




Eccentric tip-clearance distribution 



Predicted pressure time signals/spectrums at 2 locations 
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Enigma Simulations - V 
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Eccentric tip-clearance distribution 
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Blade/clearance interaction generates 10 and weak 50 unsteadiness 
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Conclusions 



Extensive analytical, computational and experimental 
program to determine the interaction of the SSME 
low-pressure fuel pump inducer with a flow liner located 
upstream 

• Two codes, Phantom and Enigma, were used to study the 
effects of operating fluid, power level, asymmetric clearance, 
duct angulation and two-phase flow on the unsteadiness 
produced from the inducer 

• The predicted results, in general, show good agreement with 
the results of water-rig and hot-fire tests 

• The predicted and experimental data show significant 
unsteadiness near several structural frequencies 

• Duct angulation and asymmetric clearance probably not 
big contributors in the crack generation process 
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